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Interdigitated  humidity  sensors  with  atomic  layer  deposited  (ALD)  coatings  of  aluminum  oxide  demon¬ 
strated  no  leakage  current  relative  to  uncoated  sensors  stored  in  the  ambient,  indicating  A1203  may  be 
used  to  limit  the  effects  of  H20  and  other  chemical  species  in  miniaturized  mechanical-  and  electronic- 
devices.  The  long  term  durability  of  such  coatings  is  not  known,  but  may  be  predicted  from  the  related 
material  characteristics.  The  modulus  and  hardness  of  A1203  were  therefore  measured  by  nanoindenta¬ 
tion  using  a  Berkovich  tip.  Because  the  coatings  are  brittle  and  possess  a  significant  tensile  stress,  the 
influence  of  film  stress  on  the  indentation  measurements  was  quantified  using  a  numerical  analysis  pro¬ 
tocol,  which  also  considered  the  effect  of  substrate  compliance.  The  film  stress  and  coefficient  of  thermal 
expansion  for  A1203  were  determined  using  the  wafer  curvature  method.  Film  stress  was  characterized 
using  thermal  cycling  up  to  500  C.  Separate  Si/Si02/Si  microcantilever  arrays  demonstrated  a  stress  vari¬ 
ation  according  to  the  thickness  of  A1203  coatings.  Fracture  toughness  was  examined  by  indentation  with 
a  cube-corner  tip;  the  estimates  are  subject  to  film  stress  and  the  material-dependent  geometry  factor. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  atomic  layer  deposition  (ALD)  technique  [1-4]  may  be 
used  to  grow  thin  metallic  or  ceramic  films.  ALD  ceramic  coat¬ 
ings  have  found  application  as  high  k  dielectrics  within  the  field 
of  integrated  circuit  technology  [5],  ALD  coatings  have  been  pro¬ 
posed  to  be  utilized  in  a  broad  range  of  applications,  including 
the  encapsulation  of  compliant  substrates  [6-8]  as  well  as  the 
surface  functionalization  of:  nano-particles  or  nano-tubes  [9,10], 
porous  films/membranes  [10,11],  and  microsystems  [12,13],  Such 
coatings  may  be  used  to  tailor  characteristics  including:  chemical 
permeation,  charge-dissipation,  surface-adhesion,  corrosion  resis¬ 
tance,  or  tribological  behavior  [13],  ALD  coatings  are  expected  to 
benefit  from  the  characteristics  unique  to  the  deposition  tech¬ 
nique:  the  resulting  films  are  continuous,  conformal,  pin-hole 
free,  and  may  be  grown  with  sub-nanometer  thickness  con¬ 
trol. 
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The  previously  unstudied  thermo-mechanical  properties  are 
essential  to  the  design  and  engineering  of  reliable  components 
containing  ALD  films.  Instrumented  indentation  [14,15]  is  a  pop¬ 
ular  technique  that  has  been  employed  to  study  a  broad  variety 
of  thin  films.  Properties,  such  as  modulus  and  hardness,  may  be 
determined  from  the  measured  load  vs.  depth  relationship  when 
a  prescribed  tip  is  impressed  into  a  specimen.  The  mechanical 
response  for  a  film  is,  however,  only  automatically  decoupled  from 
its  host  substrate  at  indentation  depths  that  are  significantly  less 
than  that  of  the  film  thickness  [16].  Owing  to  the  limitations  of 
the  technique  including:  the  tip  (capability  of  its  manufacture  and 
wear  from  use),  the  specimen  (surface  roughness,  surface  contami¬ 
nation,  and  alignment  with  respect  to  the  tip),  and  the  indent  region 
(the  evolution  if  its  initial  geometry  and  the  stress  distribution),  the 
indenter  tip  must  typically  be  pressed  into  the  specimen  by  at  least 
50  nrn  before  the  raw  modulus  and  hardness  measurements  have 
stabilized  [17].  Practically  speaking,  these  considerations  imply  the 
film  must  be  at  least  500  nm  thick  to  obtain  accurate  raw  measure¬ 
ments. 

In  addition  to  the  common  Berkovich  tip  (where  a  65.30° 
physical  angle  on  faces  results  in  an  8%  constant  strain  [15]),  a 
cube-corner  tip  may  be  used  (where  35.26°  physical  angle  results 
in  an  1 8%  applied  strain  [18]).  The  measured  E  and  H  values  as  well 
as  the  variability  may  be  affected  for  the  cube-corner  tip  for  rea- 
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sons  including  material  anisotropy  and  greater  pile-up  at  the  tip 
[19],  The  greater  applied  strain  for  a  cube-corner  tip  enhances  the 
likelihood  of  radial  crack  formation. 

The  technique  of  fracture  toughness  estimation  via  indentation 
pioneered  by  Lawn  and  others  [20-22]  requires  the  impression  site 
to  be  imaged.  The  technique  relies  on  the  assumption  that  only 
the  depth  of  the  impression  recovers  during  unloading,  whereas 
the  length  of  the  impressional  diagonals  and  the  radial  cracks 
remain  unchanged.  The  technique  originally  applied  to  shallow 
indentation  (crack  depths  <  0.4fif)  using  a  Vickers  tip  [22,23],  but 
has  been  extended  to  the  Berkovich  [24]  and  cube-corner  [25] 
tips.  Crack  depth,  which  determines  the  analysis  model,  equals 
(c-a)/2  as  corner  cracks  (generated  during  loading  or  at  the  start 
of  unloading)  typically  have  a  semicircular  profile  when  viewed 
in  cross-section  [22,23].  Eq.  (1)  relates  between  fracture  tough¬ 
ness  and  the  observed  morphology  for  deep  indentation  (crack 
depths  >0.7hf)  [26]. 

I<=  7  \  +(g[D1,D2]af(7rhf)0'5)  (1) 


In  the  equation,  here  for  system  international  (SI)  units,  I<  rep¬ 
resents  the  fracture  toughness  (MPa  m0-5);  X  the  empirical  tip 
geometry  factor  -  here  the  default  value  of  0.040  (unitless)  [25];  £ 
the  modulus  from  indentation  (Pa);  H  the  hardness  from  indenta¬ 
tion  (Pa);  it  the  mathematical  constant  (3.142);  h  the  thickness  (m); 
a  the  distance  from  the  tip  to  the  corner  within  the  impression  (m); 
c  the  distance  from  the  impression  corner  to  the  end  of  the  crack 
(m);  the  coefficient  g  accounts  for  elastic  misfit  between  the  film 
and  substrate  (unitless)  [27];  D  represents  the  Dunder’s  parame¬ 
ters  (unitless)  [28];  and  a  the  stress  (Pa).  The  subscript  f  refers  to 
the  film.  For  deep  impressions  to  qualify  for  examination,  cracks 
must  be  of  sufficient  length,  i.e.  (c-  a)  >  4hf,  which  corresponds  to 
the  condition  of  steady  state  channel-crack  propagation  [29,30].  A 
more  rigorous  model  has  recently  been  developed  for  estimating 
toughness  from  channel-cracked  films  [31],  where  a  solution  may 
be  obtained  if  the  true-  and  apparent-film  modulus,  maximum  load, 
and  other  parameters  in  Eq.  ( 1 )  are  known  for  2  or  more  indentation 
depths. 

The  curvature  of  a  film/substrate  system  is  often  used  to  exam¬ 
ine  erf  as  well  as  the  coefficient  of  thermal  expansion  (CTE).  Stoney’s 
solution  [32,33]  relates  between  curvature  and  stress  in  a  thin  film 
deposited  on  a  thick  substrate,  Eq.  (2). 


Of  = 


Aj: 

6(1  -vs)hf 


(2) 


For  curvature  vs.  temperature  profiles,  CTE  can  be  evaluated 
using  Eq.  (3)  [33], 


Ak=  6(1  -  us)hfEf(as  -  af)AT 

(1-Uf)£shs2 


(3) 


New  parameters  in  the  equations  include  k,  which  represents 
the  curvature  (m-1);  v  the  Poisson’s  ratio  (unitless);  a  the  CTE 
(ppm/°C);  and  T  the  temperature  (°C).  The  subscripts  f  and  s  refer 
to  the  film  and  substrate,  respectively.  In  practice,  the  uncoated 
substrate  is  not  perfectly  flat.  Then,  Eq.  (2)  may  be  evaluated  from 
the  difference  in  R,  the  radius  of  curvature  (m),  present  before  and 
after  deposition,  Ax  =  (l/Ra)-(l/£b)-  In  the  analysis,  the  film  is 
assumed  to  be  of  uniform  thickness,  and  rtf  is  equibiaxial  and  con¬ 
stant  throughout  hf.  If  hf  is  not  two  orders  of  magnitude  less  than  fis, 
or  for  a  stiff  film  on  a  compliant  substrate,  the  system  may  be  ana¬ 
lyzed  as  a  multilayer  composite  [34-36].  Such  analysis  is  accurate 
and  may  be  readily  performed  without  a  correction  factor  [36], 
Residual  stress  is  a  longstanding  concern  in  the  thin  film  com¬ 
munity,  and  crf  ranging  from  0.5  to  5.5  GPa  has  been  reported  for 
hf  <  500  nm  for  Cu,  Ni,  Ti,  Cr,  Mo,  Ta,  TiN,  cubic  BN,  and  Ta2N  films 


A 


OH  OH  OH 


h2o 


ch4 


Fig.  1.  Schematic  showing  sub-reaction  sequence  resulting  in  the  growth  of  AI2O3. 
Film  growth  is  realized  according  to  the  reactant  sequence  of  (a)  trimethylaluminum 
and  (b)  water. 

[37-43  ].  fff  from  3-5.5  GPa  was  directly  verified  using  X-ray  diffrac¬ 
tion  measurements  [40,42].  Stress  in  physical-  and  chemical-vapor 
deposited  films  frequently  demonstrates  a  complicated  profile 
through  the  thickness  of  the  film  [44-46],  owing  to  the  process 
of  island-formation  and  coalescence  [44-48].  Recent  study  of  the 
ALD  process  for  amorphous  materials  [49,50]  suggests  an  initial 
non-linear  growth  regime,  occurring  similar  to  island-formation, 
followed  by  steady-state  deposition  after  a  contiguous  film  has 
formed. The  mass  accumulated  duringALD  [51]  suggests  thateither 
island  formation  or  ubiquitous  growth  may  occur,  depending  on  the 
propensity  for  the  substrate  to  promote  or  hinder  film  nucleation 
(i.e.,  its  composition  as  an  oxide  or  noble  metal,  respectively). 

The  examination  here  was  motivated  by  the  use  of  ALD  to  create 
a  chemical  permeation  barrier  [8].  Regarding  that  application,  the 
long  term  durability  of  mono-  and  multi-layer  ALD  coatings  is  not 
known,  but  may  be  understood  from  the  related  material  charac¬ 
teristics.  The  properties  of  ALD  coatings  are  not  well-established, 
and  may  vary  significantly  from  those  of  bulk  material,  owing  to 
porosity  and  other  factors  [4].  The  goal  of  this  study  is  therefore 
to  examine  the  thermo-mechanical  properties  of  alumina  films, 
grown  using  the  ALD  technique.  Specifically,  the  characteristics  of 
elastic  modulus,  hardness,  film  stress,  and  CTE  are  examined  using 
indentation  and  the  wafer  curvature  method.  Drawing  on  the  lat¬ 
ter  technique,  film  stress  is  also  examined  from  the  curvature  of 
multilayer  composite  microcantilever  beams. 

2.  Experimental 

Films  were  deposited  in  a  viscous  flow  reactor  using  the  ALD 
[1-4]  technique.  Films  are  deposited  in  blanket  format,  as  the  tech¬ 
nique  does  not  require  line-of-sight  for  deposition.  The  deposition 
technique  is  based  on  a  sequence  of  two  or  more  self-limiting  reac¬ 
tions  between  vapor-phase  precursors  and  a  solid  surface.  A  simple 
recipe  for  A1203  film  growth  incorporates  the  two  half-reactions, 
(4A)  and  (4B),  where  the  asterisks  designate  the  surface  species. 

AlOH  *  +A1(CH3 )3  A10A1(CH3 )2  *  +CH4  (4A) 

A1CH3  *+H20  ->  A10H*+CH4  (4B) 

The  reactants  trimethylaluminum  (TMA,  A1(CH3)3)  and  water 
are  alternately  injected  via  nitrogen  carrier  gas.  Using  computer- 
controlled  pneumatic  valves,  the  substrate  surface  is  first  exposed 
to  TMA,  which  reacts  with  the  active  surface  sites,  Fig.  1  (a).  Then, 
after  purging  the  by-products  from  reaction  (4A),  the  surface  is 
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exposed  to  H2O.  This  reaction  regenerates  the  initial  functional 
groups,  preparing  the  surface  for  the  next  exposure  to  TMA,  Fig.  1 
(b).  The  film  is  grown  to  the  desired  thickness  by  repeating  the  AB 
sequence.  For  the  4.7 1  chamber,  the  dose  times  of  1  and  0.2  s  were 
utilized  for  TMA  and  FI2O,  respectively,  at  the  injection  pressure  of 
300mTorr.  Dosing  was  followed  by  purging  with  ultrahigh  purity 
N2  at  the  injection  pressure  of  300mTorr  for  75  s.  The  dose  and 
purge  times  of  1.5  and  120  s  were  used  specifically  for  the  leakage 
current  structures  (Fig.  4)  to  ensure  surface  exposure  and  prevent 
spurious  reaction  between  chemical  precursors,  respectively.  The 
growth  temperature  of  155°C  and  baseline  chamber  pressure  of 
650mTorr  were  used  for  all  experiments.  Direct  deposition  (with 
no  substrate  surface  treatment)  was  performed  after  a  12  h  stabi¬ 
lization  at  the  deposition  temperature. 

Indentation  of  A1203  was  performed  at  room  temperature  using 
a  commercial  instrument  (Nano  DCM,  Agilent  Technologies,  Inc.) 
equipped  with  a  diamond  Berkovich  tip.  When  measured  in  “con¬ 
tinuous  stiffness”  mode  (CSM,  Ref.  [52]),  modulus  and  hardness  can 
be  evaluated  at  discrete  instances  (unloading  events)  throughout 
the  measured  depth  range.  Once  the  tip  has  engaged  the  material’s 
surface,  the  instrument  is  capable  of  resolving  load  increments  less 
than  1  p,N,  with  displacement  resolution  less  than  1  nm.  The  instru¬ 
mented  indentation  method  is  accurate  to  within  about  5-10%  of 
the  measured  £  and  H  values  when  multiple  indentation  measure¬ 
ments  are  averaged. 

The  details  of  the  indentation  experiments  are  as  follows:  inden¬ 
tation  was  performed  to  the  depth  of  100  nm  for  the  500  nm  thick 
film  on  Si,  to  prevent  cracking  at  the  corners  as  well  as  to  limit  the 
influence  of  the  substrate.  The  test  procedure  followed  the  trape¬ 
zoidal  profile  [15]  by:  loading  at  the  constant  strain  rate  [53]  of 
0.05  (s_1 );  holding  at  the  maximum  load  for  30  s  for  stabilization: 
unloading  at  the  constant  strain  rate  of  0.05  (s_1 );  and  finally  hold¬ 
ing  at  10%  of  the  maximum  load  for  30  s  to  obtain  a  thermal  drift 
correction.  For  all  indentation  experiments,  test  locations  were  off¬ 
set  by  100  p,m  to  ensure  isolation  between  the  30  separate  sites. 

Key  details  of  the  data  reduction  of  the  raw  indentation  mea¬ 
surements  are  as  follows:  immediately  prior  to  the  tests,  the  tip 
was  calibrated  against  the  elastic  modulus  of  fused  silica  [14],  with 
the  area  coefficients  being  chosen  to  achieve  an  optimum  fit  of  the 
CSM  data  according  to  Eq.  (5). 

4 

A[/ic]  =  ^Cn(hc)2-"  =  C0h2c  +  C,hc  +  C2hl/2  +  C3hl/4  +  C4hl/8  (5) 

n= 0 


.Effective  cone 
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Fig.  2.  Schematic  summarizing  key  aspects  of  the  FEA,  including  element  types;  the 
number  of  layers;  initial  conditions;  and  boundary  conditions. 


specimen  using  Eq.  (8). 


1  1  -  v2  1  -  v2 

£eff  E{  +  Ej 


(8) 


In  the  equation,  the  subscripts  f  and  i,  refer  to  the  film  and  the 
indenter  tip,  respectively.  Et  and  Vj  were  assumed  to  be  1141  GPa 
and  0.07,  respectively.  Vf  was  assumed  to  be  0.24  [26,55]  (vs. 
Vf  =  0.13  mistakenly  identified  in  Ref.  [8]).  The  Berkovich  hardness, 
H,  was  evaluated  at  the  instant  the  tip  was  unloaded  according  to 
Eq.  (9). 


unload 


(9) 


The  raw  load,  depth,  and  harmonic  contact  stiffness  data  were 
utilized  with  a  finite  element  analysis  (FEA)  protocol  [56]  that 
inherently  incorporates  the  effect  of  the  substrate  and  allows  for 
interpretation  at  shallow  depths,  i.e„  >20  nm.  Because  the  fused 
silica  calibration  (i.e„  Eq.  (5))  does  not  ensure  a  positive  area  of 
contact  at  shallow  depths,  the  tip  was  assigned  an  ideal  spherical 
shape  for  h  <  1  nm.  For  an  assumed  elastic/perfectly  plastic  consti¬ 
tutive  behavior  profile,  the  yield  strength  and  modulus  values  are 
interpolated  from  a  fit  of  the  simulated  stiffness  and  force  relative 
to  those  measured  at  a  particular  depth,  Eqs.  (10)  and  (11). 


<7y  =  cqS  +  b\P  +  Ci 


(10) 


In  Eq.  (5),  A[h]  represents  the  depth  dependent  area  function 
(m2 )  and  C  represents  the  area  fit  coefficient(s).  The  contact  depth, 
hc,  was  determined  according  to  a  linear  fit  of  the  measured  P  vs.  h 
data,  just  after  the  tip  was  unloaded,  Eq.  (6). 

hc=h-s £  (6) 

3  unload 

New  parameters  represented  include  h,  which  represents  the 
measured  depth  ( m);  e,  the  assumed  sink-in  parameter  of  0.75  [  1 4] ; 
P,  the  applied  load  (N);  and  S,  the  “harmonic  contact  stiffness”  (i.e„ 
dP/dh)  (N/m).  The  effective  modulus  was  evaluated  at  the  instant 
the  tip  was  unloaded  according  to  Eq.  (7)  [14,15]. 


Sjri 

2/iy^M 


unload 


(7) 


New  parameters  include,  £eff,  which  represents  the  effective 
modulus  (Pa);  and  /J,  the  tip  geometry  factor  -  assumed  to  be  1.05 
based  on  Ref.  [54].  For  the  two-dimensional  axisymmetric  analy¬ 
sis,  the  effective  modulus  was  related  to  the  elastic  modulus  of  the 


E  =  a2S  +  b2P  +  c2  (11) 

£  and  ery  (therefore  H)  are  determined  iteratively  from  a  set  of 
initial  guesses  that  bound  the  converged  solution.  To  clarify,  the 
parameter  ery  represents  the  yield  strength  of  the  film  (Pa);  S,  the 
harmonic  contact  stiffness  (N/m);  P,  the  applied  load  (N);  and  £, 
the  modulus  of  the  film  (Pa).  To  facilitate  rapid  analysis,  the  mea¬ 
surements  were  binned  and  averaged  at  twenty  depths  throughout 
the  range  of  the  experiment,  where  the  fitting  coefficients  a,  b, 
and  c  are  mathematically  eliminated  at  each  depth  by  the  inter¬ 
polation  process.  In  the  protocol,  an  equibiaxial  stress  (constant 
through  the  thickness  of  the  film)  can  be  assigned  within  the  FEA 
in  order  to  examine  the  influence  of  crf  on  £  and  cry  [56],  The  values 
of£  =  168  GPa,  v  =  0.22,  and  ery  =  5270  MPa  were  determined  for  the 
substrate  (identical  to  Ref.  [57]),  from  an  analysis  of  an  indented  Si 
wafer.  The  value  of  v  =  0.24  was  assumed  for  the  AI2O3  films. 

Details  of  the  FEA,  including  the  geometry  and  boundary- 
conditions  are  summarized  in  Fig.  2.  Computation  was  performed 
using  a  commercial  code  (ABAQUS,  Dassault  Systemes  Inc.)  in 
conjunction  with  custom  front-end  utilities  [56].  The  simulated 
Berkovich  tip  (the  70.3°  equivalent  cone  geometry  for  the  two- 
dimensional  axisymmetric  condition  [15])  was  represented  using 
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a  set  of  2-node  rigid  linear-link  elements  (RAX2)  to  apply  a  dis¬ 
placement  to  an  adjoining  set  of  standard  4-node  bilinear  elements 
(CAX4R).  The  film  itself  consisted  of  a  single  layer,  10  standard 
elements  thick.  While  truncated  in  Fig.  2,  the  substrate  was  rep¬ 
resented  using  three  separate  layers  (thickness  1,  2,  and  50hf)  in 
order  to  refine  accuracy  in  the  region  near  the  tip,  while  allow¬ 
ing  the  substrate  to  approach  the  semi-infinite  condition  on  its 
far  sides.  The  modeling  of  the  substrate  as  a  semi-infinite  entity 
inherently  accounts  for  its  mechanical  compliance.  The  element 
aspect  ratio  for  the  mesh  (all  layers)  was  varied  in  accordance 
from  0.5  near  the  tip  to  5  at  the  far  boundary  in  order  to  reduce 
the  overall  computation  time.  Nodes  along  the  central-  and  far 
radial-axes  (left  and  right  of  Fig.  2)  were  radially  pinned  and 
vertically  free,  whereas  nodes  along  the  bottom  of  the  mesh  (bot¬ 
tom  of  Fig.  2)  were  radially  free  and  vertically  pinned.  An  initial 
condition  was  used  to  specify  the  film  stress.  In  the  original 
work  [56],  when  a  film  stress  is  not  present  the  4-node  linear, 
non-reflecting  one-way  infinite  element  (C1NAX4)  was  used  (bot¬ 
tom  and  right  of  Fig.  2)  to  reduce  the  overall  mesh  size,  while 
accurately  representing  the  semi-infinite  boundary  condition.  A 
contact  analysis  (assumed  friction  coefficient  of  0.2)  was  assigned 
between  the  indenter  and  film.  The  FEA  is  run  for  an  applied  dis¬ 
placement  (h  ±  1  nm  to  simulate  the  CSM  method).  The  output 
parameters  are  the  required  net  force  on  the  tip  and  the  result¬ 
ing  harmonic  contact  stiffness  of  the  tip/specimen  system.  Other 
approaches  for  the  analysis  of  the  indentation  of  film/substrate 
systems  exist,  including  solely  analytic  means  [16,58-63]  and 
numerical  analysis  [64-66],  The  instrument  measurements  are  nor¬ 
malized  relative  to  the  substrate  properties  for  the  FEA  in  Ref. 
[66]. 

A  diamond  cube-corner  tip  [15,18,19]  was  used  to  indent  the 
500  nm  thick  AI2O3  film  on  Si  up  to  400  nm,  with  a  commer¬ 
cial  instrument  (Nanoindeter  XP,  Agilent  Technologies,  Inc.)  and 
procedure  similar  to  that  described  above  (trapezoidal  loading  pro¬ 
file).  The  Nanoindeter  XP  is  capable  of  resolving  load  increments 
less  than  1  p.N,  with  displacement  resolution  less  than  1  nm.  A 
500  nm  thick  annealed  plasma  chemical  vapor  deposited  (PCVD) 
tetraethylorthosilicate  (TEOS)  SiCU  sample  [67]  was  also  indented 
for  comparison.  Here,  the  indentation  parameter  values  of  =  1.05 
and  £  =  0.75  were  also  used.  The  tip  calibration  (conducted  using 
fused  silica)  results  in  different  area  coefficients  for  the  cube-corner 
tip,  Eq.  (5).  Indentation  was  performed  adjacent  to  a  carbon  paint 
mark  so  that  the  impression  sites  could  be  located  quickly  in  a  field 
emitting  scanning  electron  microscope  (FESEM).  A  sputtered  Pt/Au 
coating  (hf  ~  1  nm)  was  added  prior  to  indentation  to  aid  FESEM 
imaging. 

AI2O3  films,  nominally  100  nm  thick,  were  deposited  on  300  p,m 
thick,  100  mm  diameter  (1  00)  Si  wafers  (University  Wafer,  Inc.)  for 
curvature  measurements  (FLX  2320-S,  Toho  Technology  Corp.).  Si 
substrates  were  utilized  because  they  are  well  characterized  and 
are  manufactured  with  high  precision.  The  use  of  thin,  large  diam¬ 
eter  wafers  compensates  the  thickness  and  modulus  of  the  AI2O3 
film,  Eq.  (2).  The  FLX  2320-S  measures  radius  of  curvature  based  on 
the  change  in  angle  (incident  vs.  reflected )  of  a  rastered  laser.  Linear 
scans  ensure  radius  measurements  to  within  2.5%  of  the  averaged 
value.  The  FLX  2320-S  has  a  temperature  controlled  chuck,  set  here 
to  record  curvature  at  each  5°C/min  increment.  To  prepare  the 
specimens  for  measurement,  A1203  was  removed  from  the  backside 
of  the  wafer  using  highly  basic  (pH~13)  solutions  of  NaOH  (from 
pellets),  hydrogen  peroxide  (30wt.%),  and  deionized  water.  The 
coating  dissolved  when  repeatedly  wiped  with  a  solution  soaked 
tissue,  as  verified  visually  from  the  indigo  appearance  of  undis¬ 
solved  A1203.  Before  and  after  wafer  curvature  measurements,  the 
film  thickness  and  index  of  optical  refraction  were  measured  a 
small  spot  spectroscopic  reflectometer  (NanoSpec,  Nanometrics, 
Inc.). 
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Fig.  3.  Micrograph  showing  the  cross-section  of  a  microcantilever  beam  coated  with 
AI2O3.  The  layers  are  labeled,  including  the  material  redeposited  underneath  by  the 
FIB. 


Similar  to  the  wafer  curvature  characterization,  AI2O3  was 
deposited  on  microcantilever  beams  in  order  to  estimate  film  stress. 
The  microcantilevers  consisted  of  laminated  polycrystalline  silicon 
(polySi),  SiC>2,  and  polySi  layers  nominally  1.5, 0.3,  and  1.0  |xm  thick, 
respectively,  Fig.  3.  The  figure  shows  an  A1203  coated  beam  after 
it  was  cut  across  its  width  using  a  focused  ion  beam  (FIB,  Nova 
Nanolab,  FEI  Company).  The  layers  present  are  labeled  according 
to  their  contrast  in  the  image,  including  the  material  redeposited 
during  the  ion  milling  process.  The  direction  of  positive  curvature 
occurs  for  tip  deflection  towards  the  substrate,  as  indicated  with 
an  arrow  in  Fig.  3.  Microcantilevers  were  fabricated  on  separate 
dice  according  to  a  standard  microsystems  technology  (SUMMiT  V, 
Sandia  National  Laboratories)  [68].  Identical  arrays  of  beams  nom¬ 
inally  20  p.m  wide  ranged  in  length  from  100  to  550  p,m  in  50  p,m 
increments.  The  microcantilevers  were  mechanically-freed  from 
sacrificial  SiC>2  layers  by  etching  in  a  solution  of  48  wt.%  hydrofluo¬ 
ric  acid  (HF)  and  Triton-X  100  surfactant  [69]  for  20  min.  As  shown 
in  Fig.  3,  the  middle  Si02  layer  of  the  beams  is  encased  in  Si,  so  that 
the  SiC>2  is  not  removed  in  HF. 

Because  the  beams  are  non-symmetric  about  their  thickness, 
they  exhibit  an  initial  curvature  of  approximately  -55  (m_1). 
In  contrast  to  the  SUMMiT  technology  [68],  the  microcantilever 
method  is  not  well  suited  for  the  MUMPs  technology  [70],  because 
the  nominal  radius  of  curvature  of  polySi/SiCU/polySi  compos¬ 
ite  structures  is  roughly  43  (m_1).  That  is,  the  tip  of  similar 
MUMPs  beams  will  contact  the  substrate  for  lengths  >200  p,m.  The 
curvature  of  the  beams  was  measured  using  an  interferometric 
microscope  (New  View  200,  Zygo  Corp.).  The  vertical  resolution 
of  the  machine  is  better  than  1  nm,  while  the  lateral  resolution  for 
thelOx  objective  at  0.75  x  magnification  is  approximately  0.89  p.m. 
The  measurement  accuracy  of  the  instrument  is  therefore  expected 
to  be  better  than  1.2%  (two  standard  deviations)  for  the  beams 
studied.  Curvature  was  measured  immediately  before  and  after 
deposition,  so  that  the  stress  in  the  coating  can  be  determined  from 
the  multilayer  composite  analysis  [34-36]. 

The  analysis  procedure  for  the  microcantilevers  is  formally 
described  in  Ref.  [36];  the  final  step  is  to  relate  between  micro¬ 
cantilever  curvature  and  stress,  Eq.  (12).  New  parameters  in  the 
equation  include  A,  the  constant  coupling  axial  extension  between 
the  layers  (N);  D,  the  bending  coupling  constant  (N  m2);  B,  the  con¬ 
stant  coupling  between  extension  and  bending  (N  m);  N,  the  term 
for  the  laminate  force  (N);  and  M,  the  term  for  the  laminate  bending 
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Fig.  5.  Raw  indentation  results  for  500  nm  thick  AI2O3  film  on  Si. 


Fig.  4.  Leakage  current  measurements  for  an  interdigitated  leakage  current  sensor. 
The  inset  contains  an  optical  micrograph,  showing  the  test  structure  from  the  top. 
The  applied  electrical  conditions  (//V)  indicated  in  the  inset. 


moment  (N  m).  The  coefficients  A ,  B,  and  D  vary  with  the  E,  v,  and  h 
of  each  of  the  component  layers.  The  terms  N,  and  M  vary  with  the 
E,  v ,  h,  a,  and  a  of  the  component  layers  as  well  as  AT.  The  analy¬ 
sis  was  simplified  to  the  one-dimensional  condition  (layer  width  is 
eliminated)  because  the  beam  length  was  significantly  greater  than 
the  width. 


(- BN)  +  {AM ) 

(AD)-m 


(12) 


3.  Results 


3.1.  Leakage  current  assessment 

ALD  surface  coatings  are  of  interest  as  a  moisture  barrier 
enabling  wafer  level  encapsulation  (WLE)  of  microelectromechani¬ 
cal  systems  (MEMS).  The  WLE  is  formed  by  first  globally  depositing 
a  sacrificial  (SAC)  layer  and  then  locally  patterning  the  SAC  layer 
at  the  MEMS  devices.  Then  an  encapsulation  structural  layer,  typ¬ 
ically  an  oxide  (LOX),  is  deposited  over  the  sacrificial  layer.  The 
LOX  is  patterned  and  etched  to  expose  the  SAC  layer  that  is  subse¬ 
quently  removed  by  a  sacrificial  release  process.  In  the  next  steps, 
the  release  etch  holes  are  plugged  and  a  sealing  barrier  layer  is 
added.  In  the  experiments  here,  supplemental  ALD  barrier  films 
were  added  for  evaluation  in  the  absence  of  the  final  sealing  barrier. 

Evaluation  of  the  WLE  is  performed  using  a  leakage  current 
monitor.  The  leakage  current  monitor  is  comprised  of  a  series  of 
interdigitated  conductive  fingers,  as  shown  in  Fig.  4  (inset).  The 
interdigitated  structure  has  traditionally  enabled  dew  point  detec¬ 
tion  during  plastic  package  qualification  [71-73].  After  calibration 
in  a  temperature-  and  moisture-controlled  environment,  the  same 
structure  may  function  as  a  humidity  sensor,  whether  encapsulated 
by  a  permeable  polymeric  layer  [74]  or  direct  exposed  to  the  envi¬ 
ronment  [75,76],  In  the  direct  ambient,  monolayers  of  water  at  the 
surface  may  facilitate  electric  current. 

The  leakage  current  monitor,  Fig.  4,  is  fabricated  at  a  commer¬ 
cial  foundry  from  an  Al  layer  with  a  line  and  space  of  5  p,m.  The 
leakage  current  monitor  is  encapsulated  at  the  wafer  level  so  that 
ALD  surface  coatings  can  be  added  directly  on  the  WLE  to  act  as  a 
barrier  to  moisture.  This  approach  to  validate  the  WLE  has  been  pre¬ 
viously  used  to  examine  thin  film  [74,77]  and  wafer  bonded  [75,76] 
encapsulation  schemes. 

As  shown  in  Fig.  4,  current/voltage  traces  were  recorded  in 
the  ambient  environment  for  uncoated  sensors  using  a  semicon¬ 
ductor  parameter  analyzer  (Agilent  Technologies,  HP  4155B)  with 


worst  case  accuracy  of  ±0.73%.  In  contrast  to  the  uncoated  (“ambi¬ 
ent”)  condition,  no  signal  was  detected  within  the  resolution  of  the 
instrument  (1  x  10-14A)  for  the  same  sensors  immediately  after 
they  were  baked  at  200  °C  for  24  h.  No  current  was  detected  for  sen¬ 
sors  coated  with  25  nm  A1203, 50  nm  A1203,  or  25  nm  A1203  +  25  nm 
Al203/Zn0.  A  slight  leakage  current,  however,  was  detected  for 
sensors  coated  with  25  nm  A1203  +  60nm  Si02. 

Fig.  4  demonstrates  the  utility  of  ALD  coatings  for  encapsulation 
in  applications  where  a  hermetic  package  may  not  be  physically 
possible  or  cost-effective.  A  monolayer  of  A1203  may  provide  mois¬ 
ture  protection,  preventing  electrical  shorting  or  spurious  corrosion 
currents.  The  h f  of  25,  and  50  nm  were  explored  here  for  A1203 
to  ensure  adequate  breakdown  voltage  [78].  While  A1203  is  prone 
to  corrode  in  strongly  acidic  or  basic  environments  [79],  A1203 
may  be  used  as  the  base  layer  for  other  inert  materials  -  such 
as  chemical  vapor  deposited  Si02  [79].  Alternately,  an  electrically 
conductive  Al203/Zn0  composite  may  be  deposited  onto  A1203  for 
the  purpose  of  electric  charge  dissipation.  For  Al203/Zn0,  the  top¬ 
most  layer  is  a  ceramic  composite  [80]  grown  by  alternating  the 
deposition  chemistry  (8  cycles  of  A1203  and  5  cycles  of  ZnO),  result¬ 
ing  in  a  50%  volumetric  mixture.  ZnO  was  grown  here  using  the 
diethylzinc  (DEZ,  Zn(CH2CH3  )2 )  and  H20  chemistry  system  [80].  All 
of  the  aforementioned  coatings  proved  very  successful  in  reducing 
the  leakage  current  relative  to  the  uncoated  configuration  in  the 
ambient  condition.  The  distinct  profile  for  the  25  nm  A1203  +  60  nm 
Si02  coating  may  be  caused  by  residual  stress  related  cracking;  the 
mechanism  enabling  leakage  current  was  not  specifically  investi¬ 
gated. 

3.2.  Instrumented  indentation 

To  understand  the  durability  of  the  coatings  in  Fig.  4,  the  key 
thermo-mechanical  properties  were  investigated.  Binned  and  aver¬ 
aged  indentation  data  is  shown  in  Fig.  5,  including  the  variation  in 
P/S2  and  load  with  indentation  depth.  P/S2  [81]  is  seen  to  asymp¬ 
totically  converge  at  about  h  =  50  nm.  The  initial  P/S2  profile  (from 
fi  =  0-40nm)  identifies  that  the  depth  of  contact  is  not  sufficient 
to  be  readily  represented  using  the  standard  analysis  procedure, 
e.g.,  because  of  the  imperfect  geometry  of  the  tip  [14,82],  The  more 
uniform  behavior  for  h  >  40  nm  results  from  the  mechanical  char¬ 
acteristics  of  film  followed  by  the  film/substrate  system.  The  slope 
of  the  P/S2  profile  at  greater  depths  (e.g.,  h  =  100  nm)  may  be  used 
to  assess  the  contribution  of  the  substrate  relative  to  that  of  the 
film.  While  the  specimen  was  further  indented  up  to  /tf/5,  the  val¬ 
ues  of  £=183.9 ±6.46  and  H  =  1 1.6± 0.7 GPa  were  determined  for 
/if  =  50  ±  5  nm.  The  P  vs.  h  profile  in  Fig.  5  is  typical  of  a  ceramic  mate¬ 
rial  [14,15]  and  is  comparable  to  that  previously  observed  for  ALD 
A1203  [83,84],  Because  the  raw  £f  measurement  slightly  exceeds 
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Fig.  6.  Indentation  results,  analyzed  for  the  500  nm  thick  A1203  coating  on  Si. 


£s,  and  the  film  stress  of 474  MPa  has  been  previously  measured  for 
ALD  A1203  [83  ],  the  indentation  data  was  subjected  to  the  numerical 
analysis  protocol  [56]. 

As  shown  in  Fig.  6,  the  indentation  data  was  evaluated  over  a 
broad  range  of  trf  conditions.  To  clarify,  E  and  H  in  Fig.  6  were  eval¬ 
uated  from  the  indentation  data  in  Fig.  5  as  a  function  of  the  Of 
that  might  be  present.  As  described  later,  a f  was  independently 
verified  to  vary  with  temperature.  In  the  figure,  E  is  seen  to  vary 
monotonically.  H  further  depends  on  inelastic  material  behavior 
and  varies  more  complexly  with  stress  in  Fig.  6.  Results  for  af  =  0  and 
<7f  =  500  MPa  are  given  for  reference  in  Table  1 .  The  yield  strength 
and  corresponding  strain,  ey,  are  provided  in  the  table  for  the 
assumed  linear  elastic  limit. 

To  begin  discussion  of  the  indentation  results,  Fig.  6  and  Table  1 
identify  minor  variation  with  a f  over  the  range  of  500  MPa.  For 
AI2O3,  a  3.8  and  2.1%  variation  with  erf  is  observed  for  E  and  H, 
respectively.  Regarding  variation  in  H  with  Cf,  the  <Tf  of  500  MPa 
comprises  10.2%  of  cry,  Table  1,  whereas  the  crf  of  2000  MPa 
approaches  41%  of  cry.  During  applied  tensile  stress,  the  film 
becomes  more  inclined  to  flow  during  indentation,  decreasing  its 
hardness.  The  geometry  of  the  residual  impression  is  accurately 
rendered  in  FEA.  The  influence  of  geometry  on  E  (occurring  through 
S  as  well  as  the  corresponding  area  of  contact  in  Eq.  (7))  is  therefore 
considered  in  Fig.  6  and  Table  1. 

Regarding  the  E  and  H  values  in  Table  1,  they  are  similar  to 
those  previously  measured  for  thin  A1203  films.  A  survey  of  the  lit¬ 
erature  in  Ref.  [83]  identifies  100<£<272GPa  and  8<H<9.6GPa 
for  evaporated,  chemical  deposited,  and  vapor  deposited  A1203. 
A  recent  study  specific  to  ALD  A1203  identified  the  E  and  H  val¬ 
ues  of  220  and  10.5  GPa,  respectively,  using  an  analytic  model  to 
account  for  the  film/substrate  geometry  [84].  In  comparison,  the 
Hill  modulus  for  untextured  polycrystalline  a-AI203  (“corundum”) 
is  402.7  GPa  [85],  while  the  measured  E  for  bulk  amorphous  A1203 
is  372  GPa  [86].  The  modulus  for  ALD  A1203  is  decreased  relative 
to  the  bulk  value  because  of  its  amorphous  nature  [4[.  That  is, 
the  density  of  ALD  A1203  is  3.0g/cm3  [87],  whereas  that  of  bulk 
monocrystalline  a-Al203  is3.96g/cm3  [88].  While  the  complicated 
three-dimensional  strain  field  invoked  during  indentation  does 
not  readily  distinguish  anisotropic  property  variation  [89,90],  the 
isotropic  structure  of  ALD  A1203  [4]  prohibits  such  variation  here. 


Table  1 

Summary  of  analysis  results  (average  ±  2  SD)  for  500  nm  thick  A1203  film  on  Si. 


Of  (MPa) 

Ef  (GPa) 

Hf  (GPa) 

oy  (MPa) 

£y(%) 

0 

195.3  ±  20.1 

6.65  ±  0.80 

4,882  ±  679 

2.5 

500 

202.7  ±  21.8 

6.79  ±  0.79 

5,177  ±  644 

2.6 

500  nm 


300  nm 


Fig.  7.  FESEM  images  of  the  residual  impression  remaining  after  indentation  to 
h  =  400  nm  using  a  cube-corner  tip  for  500  nm  thick  (a)  A1203,  and  (b)  Si02  coatings 
on  Si. 

The  numerical  analysis  protocol  also  importantly  identifies  the 
influence  of  the  substrate,  which  was  not  previously  considered  in 
Refs.  [12,83].  The  mechanical  compliance  of  the  substrate  (occur¬ 
ring  because  Es  <£f)  increases  £f  by  11.4  MPa  (6.2%)  over  the  raw 
measured  value.  In  comparison,  an  additional  increase  of  7.4  MPa 
for  the  corrected  £f  is  identified  for  <rf  =  0  ->  500  MPa.  Similar  trends 
are  observed  for  H f,  although  the  base  value  of  6.65  GPa  (which 
accounts  for  substrate  compliance)  is  less  than  that  in  previous 
studies  [12,83,84].  While  the  P/S2  vs.  h  profile  in  Fig.  5  suggests 
that  the  film  is  largely  decoupled  from  substrate,  the  minor  effect 
of  substrate  compliance  might  have  been  anticipated  from  the  slope 
of  a  least  squares  fit,  applied  in  the  vicinity  of  /if/10.  In  comparison 
to  the  raw  measured  values,  the  FEA  rendered  valid  (converged) 
data  at  30  nm  (0.06  x  fif). 

3.3.  Indentation  toughness 

Representative  examples  of  the  residual  impression  remaining 
after  cube-corner  indentation  are  shown  in  Fig.  7  for  the  separate 
ALD  A1203  and  PCVD  Si02  coatings.  Radial  cracking  was  observed 
at  the  corners  of  the  A1203  impressions,  beginning  at  indentation 
depth  of  100  nm  (for  data  sets  of  h  =  50, 100,  200,  300,  and  400  nm). 
Cracks  for  the  Si02  impressions  never  became  of  sufficient  length  to 
be  considered  for  analysis,  i.e.,  c/a  >2.5  for  h<0.4hf  or  (c-  a)>  41if 
for  h<0.4/if.  Only  the  fi=400nm  data  for  A1203  became  of  suffi¬ 
cient  length  to  be  analyzed  as  a  channel-cracked  film.  In  no  cases 
did  the  film/substrate  systems  spall,  such  that  would  invalidate  the 
analysis.  Analysis  was  conducted  from  the  images  of  the  impres¬ 
sion  site,  as  there  was  no  indication  of  the  onset  of  cracking  in 
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Fig.  8.  Stress  results  over  the  first  3  consecutive  thermal  cycles  for  100  nm  AI2O3 
deposited  on  Si.  The  specimen  geometry  for  the  wafer  curvature  measurements  is 
shown  in  the  inset. 
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the  indentation  data  profiles  according  to  the  criteria  of  P  vs.  h 
[18];  load/stiffness  squared  (P/S2)  vs.  h  [18,81];  Kexp  vs.  h  [18,91]; 
or  reduced  stiffness  (Sc)  vs.  h  [18],  For  A  =  0.040,  £f  =  195.3 GPa, 
vf  =  0.24,  Hf  =  6.65  GPa,  hf  =  500  nm,  £s  =  1 61 .8  GPa,  and  vs  =  0.22,  the 
nominal  Kjc  of  0.76  ±  0.18  MPa  m0-5  was  determined  for  <rf  =  0.  If  Of 
was  allowed  to  vary,  the/Cic  of  1.89  MPa  m0-5  (measured  previously 
in  Ref.  [8])  resulted  for  the  tensile  stress  of  838  MPa. 

A  recent  review  discussing  the  limitations  of  estimating  tough¬ 
ness  via  indentation  [92]  reminds  that  the  legacy  of  the  method 
owes  to  empirical  means  and  not  a  theoretical  origin.  In  particu¬ 
lar,  X  may  vary  according  to  the  specimen  material.  Comparison 
against  the  Si02  specimen  also  suggests  a  significant  tensile  a( 
for  the  AI2O3  film,  such  that  would  aid  corner  cracking.  (The 
Si02  specimen  has  the  same  geometry  but  different  composi¬ 
tion.)  To  explain,  the  Si02  [67]  was  annealed  at  high  temperature, 
which  would  relax  Of,  making  it  less  prone  to  cracking.  Eq.  (2) 
and  (3)  may  be  combined  to  estimate  crf  for  the  Si02  specimen, 
i.e„  <Tf  =  (fif/(l  -  Vf))(af- as)AT.  The  values  of£f  =  70GPa,  Vf  =  0.18, 
o?f  =  0.4ppm/°C,  and  cfs  =3.0ppm/°C  imply  crf  =  -11 7  MPa  for  Si02. 
Compressive  crf  (invoked  after  the  Si02  cooled  from  the  glassy  state 
(Tg~  550  °C)  to  the  ambient-temperature)  would  hinder  crack  for¬ 
mation.  While  the  value  of  X  =  0.040  was  utilized  for  the  analysis 
above,  the  appropriate  value  is  not  well  established  in  the  literature 
for  the  cube-corner  tip  geometry.  Eq.  ( 1 )  renders  X  =  0.068  ±  0.002 
for  /< ic  =  1 .89  MPa  m0-5  if  the  tensile  crf  of  422-470  MPa  is  assumed 
for  the  A12C>3  film  (see  Fig.  8,  below).  Further  study  is  warranted, 
because  the  technique  in  Fig.  7  is  limited  by  the  material  dependent 
X  and  process  specific  crf  values. 


3.4.  Wafe rcurvatu re 

Fig.  8  shows  the  stress  present  in  A12C>3,  measured  using  the 
wafer  curvature  method  for  three  consecutive  thermal  cycles. 
The  specimen  was  specifically  heated  to  275  °C  and  then  actively 
cooled  to  25  °C.  The  separate  thermo-mechanical  profiles  cannot 
be  readily  distinguished  in  Fig.  8.  The  initial  crf  of  422  ±  21  MPa 
(present  at  the  beginning  of  the  first  thermal  cycle)  was  deter¬ 
mined  by  comparing  coated  and  uncoated  wafers.  During  the 
all  three  thermal  cycles,  stress  varied  linearly  with  temperature. 
Similar  response  was  observed  in  a  separate  specimen  heated 
to  500  °C.  The  radius  of  curvature  ranged  from  -150  to  -340  m 
during  the  characterization  in  Fig.  8.  Eq.  (3)  (applied  to  the  least- 
squares  fit  in  Fig.  8)  identifies  the  CTE  of  4.2  ±0.1  ppm/°C  for 
ALD  A1203.  A  2.5-nm  thick  native  oxide  [93]  is  expected  to  be 


present  at  the  interface  between  Al203  and  Si.  This  interlayer  is 
expected  to  have  negligible  influence  here,  because  it  is  symmetric 
(present  on  both  the  top-  and  bottom-Si  surfaces)  and  is  signif¬ 
icantly  thinner  than  /if.  The  thickness  of  A12C>3  was  determined 
by  the  reflectometer  to  be  101.9  and  102.3  nm  before  and  after 
three  thermal  cycles,  respectively.  Similarly,  the  refractive  indices 
of  1.641  and  1.640  were  measured  for  A12C>3  before  and  after  the 
experiment. 

The  CTE  of  4.2  ppm/°C  is  less  than  that  recently  measured  for 
amorphous  A1203  thin  films,  where  a  =  7.1  [26]  and  5.0  ppm/°C  [55], 
For  bulk  Al203,  the  CTE  of  5.1  and  6.6ppm/°C  has  been  identified 
for  crystalline  [94]  and  amorphous  [86]  material.  The  lesser  cq  for 
ALD  A1203  may  owe  to  its  density,  which  is  roughly  3/4  that  of  fully 
dense  crystalline  material  [87,88]. 

The  estimates  of  CTE  from  Fig.  8  allow  the  contribution  of  ther¬ 
mal  stress  to  be  evaluated.  The  values  of  £f  =  195.3  GPa,  vf  =  0.24, 
af  =  4.2  ppm/°C,  and  as  =3.0ppm/°C  imply  the  tensile  stress  of 
35  MPa  is  generated  by  cooling  from  the  deposition  temperature 
of  155  °C  to  the  ambient.  The  contribution  of  thermal  misfit  is 
therefore  estimated  to  be  5-10%  of  Of  =  422  MPa.  The  remaining 
stress  is  attributed  to  other  factors,  such  as  film  growth.  For  A1203, 
the  parameter  of  temperature  (but  not  pressure)  affects  the  rate 
of  growth  according  to  the  population  of  active  surface  sites  [3], 
Deposition  temperature  might  therefore  be  used  to  manipulate  af. 
The  added  tensile  stress  of  67  MPa  is  predicted  for  the  temperature 
change  from  275  to  25  °C  in  Fig.  8. 

The  three  overlapping  k  vs.  T  profiles  in  Fig.  8  suggest  ALD  A12C>3 
is  stable  over  the  temperature  and  time  scale  of  the  experiment. 
To  clarify,  material  stabilization  occurring  during  thermal  cycling 
is  common  in  thin  films.  Stabilization  is  motivated  by  mechanisms 
including  grain  growth,  recrystallization,  or  (inter-/intra-)  diffusion 
[37,95,96].  Material  stabilization  typically  results  in  a  complicated 
k  vs.  T  profile  during  the  first  thermal  cycle,  k  is  then  linear  with  T, 
so  long  as  the  previous  T-  or  t-conditions  are  not  exceeded.  The  hf 
and  11  values  before  and  after  thermal  cycling  were  identical  (within 
the  instrument  accuracy),  further  suggesting  that  A12C>3  is  stable. 
The  curvature  of  microcantilever  beams  coated  with  A1203  (as  in 
Fig.  3)  was  found  in  Ref.  [97]  to  be  unchanged  over  4  months  in  the 
ambient  environment,  additionally  suggesting  that  A1203  is  stable 
with  time. 

3.5.  Microcantilever  curvature 

Separate  Al203  coatings  of  different  thickness  were  grown  on 
separate  arrays  of  polySi/Si02/polySi  microcantilever  beams.  For 
the  microcantilevers  of  known  thickness  (Fig.  3),  the  parameters  of 
£f=  195.3  GPa,  vf  =  0.24,  af  =  4.2  ppm/°C,  £si  =  161.8  GPa,  vsi  =  0.22, 
asi  =  3.0ppm/°C,  Esic,2  =  70.0  GPa,  vsic,2=0.18,  asio2  = 

0.4  ppm/°C,  and  AT=130°C  were  used  in  Eq.  (12)  to  esti¬ 
mate  Of.  The  residual  stress  in  the  two  polySi  layers  in  Fig.  3 
was  determined  to  be  -2.8  and  -4.3  MPa,  respectively,  using  the 
pointer  structure  and  procedure  described  in  Ref.  [98].  This  may 
include  the  thermal  contraction  of  polySi  occurring  between  the 
deposition  and  ambient  conditions.  Residual  stress  in  the  Si02 
layer  was  assumed  to  be  0.  The  erf  values  in  Table  2  were  deter¬ 
mined  from  the  difference  in  curvature  before  and  after  coating. 
The  critical  tensile  stress  required  for  channel-crack  propagation, 
erc.  is  also  shown  in  Table  2.  <rc  was  estimated  for  the  Al203/Si 
system  according  to  Beuth’s  model  [27]  using  the  parameters  of 
£f=  195.3  GPa,  vf  =  0.24,  £s  =  161.8GPa,  vs  =  0.22  and  the  mode  1 
fracture  toughness  of  1.89  MPa  m0-5  forAl203  [8[.  The  erf  estimates 
in  Table  2  never  exceed  ac.  Thermal  misfit  at  the  Al203/polySi 
interfaces  is  inherently  incorporated  into  the  results  in  Table  2, 
based  on  the  temperature  change  between  the  deposition  and 
ambient  conditions.  The  contribution  of  thermal  misfit  became 
increasingly  significant  as  lif  was  decreased,  affecting  crf  by  25%  for 
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Table  2 

Estimate  of  residual  in-plane  stress  for  A1203,  from  comparison  of  curvature  before 
and  after  coating  microcantilever  beams.  The  implied  through-thickness  gradient 
may  come  from  the  geometry  of  the  component  layers  as  well  as  the  morphology 
of  the  Al203/Si  interfaces. 


hf,  ALD  AI2O3  (nm) 

Measured  A K  (m-1 ) 

o-f  (MPa) 

<tc  (MPa) 

10 

2.18  ±  1.04 

10140  ±  1120 

13140 

25 

4.25  ±1.11 

4980  ±  500 

8310 

50 

9.51  ±  2.60 

3770  ±  600 

5880 

100 

10.51  ±  1.36 

2200  ±  180 

4160 

125 

11.98  ±  1.56 

2010  ±  180 

370 

500 

4.07  ±  0.51 

470  ±  50 

1860 

/if=10nm.  In  Fig.  8  and  Table  2,  erf  is  importantly  assumed  to  be 
constant  through  hf. 

The  curvature  values  in  Table  2  are  troubling,  not  only  because 
they  imply  the  analysis  of  Fig.  8  is  limited,  but  also  because  they 
are  inconsistent  with  the  measurements  in  Ref.  [83],  The  a f  val¬ 
ues  in  Table  2,  however,  should  not  be  interpreted  literally.  The 
gradients  implied  in  the  table  likely  owe  more  to  the  mechan¬ 
ical  structure  of  the  microcantilevers  than  the  characteristics  of 
the  component  materials.  To  explain,  the  curvature  of  the  micro¬ 
cantilevers  occurs  as  the  net  deformation  relieving  the  different 
stresses  within  each  of  the  different  thin  films.  A  through-thickness 
strain  gradient  within  the  A1203  layer  itself  (incurred  during  the 
process  of  the  formation  of  the  AI2O3  layer)  is  not  required  to  render 
the  results  in  Table  2. 

For  the  four  films  in  Fig.  3,  crf  may  occur  because  of  stress 
remaining  from  the  process  of  film  formation;  the  morphology  of 
the  different  film  surfaces;  lattice  mismatch  at  the  interfaces;  and 
CTE-misfit  at  the  interfaces.  Significant  <jf  may  develop  in  chemical 
vapor  deposited  (CVD)  films  within  the  first  10-30  nm  of  material, 
owing  to  the  process  of  island  formation  and  coalescence  [44-46]. 
Because  the  native  oxide  on  Si  will  readily  facilitate  the  initial  gen¬ 
eration  of  hydroxyl  groups  at  the  surface  (Eq.  (4A)),  the  island 
formation  mechanism  is  not  expected  for  the  ALD  A^Oa/polySi 
system. 

The  morphology  of  the  different  polySi  film  surfaces  may  greatly 
contribute  to  the  deformation  of  the  microcantilevers.  Crevices 
are  deeper  (120  vs.  55  nm)  and  the  surface  roughness  is  greater 
(15.6±  12.8  nm  vs.  5.4±4.0nm,  forRavg  ±krms)  at  the  top  and  bot¬ 
tom  surfaces  of  polySi,  respectively  [99].  The  different  porosity  at 
the  polySi  surfaces  dictates  that  a  different  volume  of  AI2O3  will 
exist  adjacent  to  the  top  and  bottom  A^C^/polySi  interfaces.  The 
different  geometry  present  at  the  interfaces  therefore  allows  a  CTE 
motivated  bending  moment.  In  contrast,  the  surface  of  polished  Si 
wafers  lacks  grain  boundary  grooves.  Further,  the  roughness  of  Si 
wafers  is  typically  on  the  order  of  0.3  ±  0.4  nm  [69],  i.e.,  comparable 
to  the  ALD  growth  increment.  While  the  argument  of  a  mechanical 
bending  moment  enabled  by  the  surface  morphology  in  addition  to 
the  geometry  and  location  of  the  component  layers  is  not  the  only 
possible  explanation,  it  does  allow  for  significant  cr{  in  the  absence 
of  a  through-thickness  strain  gradient  for  an  ALD  Al2C>3/Si  wafer 
system. 

4.  Conclusions 

Interdigitated  leakage  current  sensors  with  atomic  layer 
deposited  (ALD)  coatings  of  (a)  AI2O3,  (b)  Al203/Si02  layers,  or 
(c)  Al203/Zn0  layers  demonstrated  no  leakage  current  relative 
to  uncoated  sensors  stored  in  the  ambient.  This  identifies  that 
ALD  AI2O3  may  be  used  to  limit  the  effect  of  H2O  and  other 
chemical  species.  To  enable  durability  predictions,  the  thermo¬ 
mechanical  properties  of  ALD  AI2O3  were  therefore  investigated 
using  indentation,  wafer  curvature,  and  microcantilever  curvature 
measurements.  Key  results  include: 


For  zero  film  stress,  the  elastic  modulus  of  AI2O3  was  deter¬ 
mined  to  be  195.3  GPa  using  indentation,  while  the  Berkovich 
hardness  was  6.65  GPa.  A  2-4%  variation  in  modulus  and/or  hard¬ 
ness  according  to  film  stress  was  identified  using  a  numerical 
analysis  protocol.  The  variation  with  stress  is  about  half  of  that 
associated  with  the  mechanical  compliance  Si  substrate,  rendering 
the  final  corrected  modulus  and  hardness  of  202.7  and  6.79  GPa, 
respectively. 

The  coefficient  of  thermal  expansion  for  AI2O3  was  determined 
to  be  4.2  ppm/°C  using  the  wafer  curvature  method.  The  total  stress 
present  after  deposition  was  422  MPa.  The  tensile  thermal  stress 
invoked  during  cooling  from  the  deposition  temperature  to  the 
ambient  was  identified  to  be  35  MPa,  i.e.  5-10%  of  the  film  stress. 
Stress  decreased  linearly  during  thermo-mechanical  loading,  but 
was  not  changed  in  the  ambient  condition  despite  repeated  ther¬ 
mal  cycling  or  temperature  excursion  to  500  C.  The  thickness  and 
index  of  refraction  were  unchanged  after  thermal  cycling,  suggest¬ 
ing  A1203  is  stable  throughout  the  temperature  range  examined. 

Separate  Si/Si02/Si  microcantilever  arrays  demonstrated  a  cur¬ 
vature  variation  according  to  the  thickness  of  AI2O3  coatings.  The 
shape  change  is  expected  to  depend  on  the  geometry  and  location 
of  the  component  layers  as  well  as  the  different  porosities  at  the 
top-  and  bottom-Si  surfaces. 

For  the  film  stress  of  422  MPa  and  mode  1  fracture  toughness  of 
1.89  MPa  m0-5,  the  material-dependent  geometry  factor  for  cube- 
corner  indentation  may  be  as  great  as  0.068.  In  the  absence  of 
previously  determined  stress  or  toughness  characteristics,  the  raw 
toughness  of  0.74  MPa  m0-5  or  tensile  film  stress  of  838  MPa  were 
independently  identified  here. 
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